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Pade approximants to truncated post-Newtonian neutron star models are constructed. The Pade 
models converge faster to the general relativistic (GR) solution than the truncated post-Newtonian 
ones. The evolution of initial data using the Pade models approximates better the evolution of full 
GR initial data than the truncated Taylor models. In the absence of full GR initial data (e.g., for 
neutron star binaries or black hole binary systems), Pade initial data could be a better option than 
f*"^ ' the straightforward truncated post-Newtonian (Taylor) initial data. 
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I. INTRODUCTION 

Compact binary systems of neutron stars and black holes inspiralling under gravitational radiation reaction are 
one of the most promising sources of gravitational waves for kilometre arm length interferometric gravitational wave 
detectors like the LIGO and VIRGO. The inspiral phase is best described by a post-Newtonian approximation which 
should eventually break down and the final merger and coalescence may be only accessible via numerical integration 
■^j- ■ of Einstein's equations. A major obstacle in the numerical studies of such systems is the non-availability of physically 
Q\ \ satisfactory initial data. Given the constraints on computational resources, one would like to start the numerical 
integration as close as possible to the final coalescence phase, using the initial data obtained by matching on to the 
known analytical results of inspiral. One of the suggestions in this direction is to use the analytical post-Newtonian 
results of inspiral to provide initial data (e.g. jjj]- 0]) for the numerical integration of the fully general relativistic 
system. 

Before attempting to apply the above strategy to the complicated compact binary system, as a preliminary step, 
Shinkai Q| has constructed a single neutron star model using the post-Newtonian approach and concluded that the 
truncated second order post-Newtonian approximation is close enough to describe a general relativistic single star. 
£~ I ■ The truncated post-Newtonian series used above is essentially a Taylor expansion in the three small parameters 
A = p/(p t c 2 ), B = 4npr 3 / (mc 2 ) and C — 2Gm/ '(rc 2 ) where p is the pressure, pt, the total energy density of the 
system, r the radial coordinate, and m(r) the mass contained in a sphere of radial coordinate r. We shall refer to 
the post-Newtonian truncated models alternatively as Taylor models. Recently, in a related context of gravitational 
wave phasing, the slow convergence of straightforward Taylor approximants has been critically investigated [|| ||. It 
was shown that new approximants, with much improved convergence properties, may be constructed for gravitational 
wave data analysis applications using, as an important tool, Pade techniques to estimate the relevant functions from 
only the first few terms in their perturbative post-Newtonian expansion. This approach has been systematically 
extended in a series of publications [0 and most recently has been used to go beyond the adiabatic approximation 
to inspiral and provide an analysis of the transition from the inspiral to the plunge in binary black hole coalescences 
Q . This work also provides for the first time initial dynamical data (positions and momenta) for binary black holes 
starting to plunge, so that there is less than an orbit left to evolve. In view of this experience, in this report, we 
construct Pade approximants to the truncated post-Newtonian (Taylor) neutron star models discussed by Shinkai || 
and investigate their performance. We show that the Pade models are better than the truncated post-Newtonian 
models and converge faster to the exact general relativistic solution. Further, we also show that the evolution of a 
general relativistic (single) star is described much better by Pade initial data as compared to the truncated Taylor 
initial data of the same order. 

In the next section we discuss the TOV equation and the truncated post-Newtonian models. In Section III we 
introduce the standard Pade approximant, discuss its applicability to the present problem, and adapt it to construct 
an appropriate generalisation for the case of 3-small parameters. Pade models corresponding to the 'Taylor' truncated 
models are constructed. In the last section we discuss the results and summarise our conclusions. The Appendix lists 
the more involved formulas for the 3PN Pade models. 
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II. TOV EQUATION AND TRUNCATED POST-NEWTONIAN (TAYLOR) NEUTRON STAR MODELS 



In general relativity the metric of a spherically symmetric static star can be written as 

ds * = - e Mr) dt 2 + e 2A(r) dr 2 + f2( ^2 + r 2 ^2 (2J) 

and the equations of hydrostatic equilibrium - the Tolman-Oppenheimer-Volkoff (TOV) equations - obtained (in 
gcometrised units: G = c = Mq = 1) from the Einstein field equations, for a given fluid distribution specified by an 
adiabatic equation of state p = p(pt), are given by: 

dm o ,„ „, 

4tr 2 ft , (2.2) 



d P _ mp t f 1+ p\ LVU^toV^ (23) 



dr 

dr r 2 \_ Pt J V m 

dv m ( Airpr 3 \ ( 2m\ 

' 1 + — — 1 • 2.4 

m I \ r 1 



Ira 



dr r 2 

m is given in terms of metric component A as 

e2A =^-vJ • (2 - 5) 

The above set of equations are integrated from the center (r = 0) to the boundary of the star, with the initial 
conditions m(r = 0) = 0;p t (r = 0) = p tc and v(r = 0) = v c . v c is rescaled appropriately such that it matches with 
the exterior Schwarzschild solution at the boundary. The radius of the star, R, is characterised as the radius where 
density p t {r — R) drops to 10 6 gm/cc (approximately 0(1O~ 10 ) in geometrised units). 



Schematically, equations (2.3) and (UA) can be written as 



d f = ^ [l + A){l + B)(l-C)~\ (2.6) 
dr r z 

*L = ™(i + B )(l C)~\ (2.7) 
dr r z 

Assuming that A, B and C are of comparable orders of smallness and making a Taylor series expansion in A, B, and 
C around the origin, the above equations to third post-Newtonian order then yield, 

|? = -^(1 + A + B + C + AB + BC + CA + C\ 

1PN 2PN 

+ ABC + AC 2 +BC 2 + C 3 +0(4)---), (2.8) 



3PN 

dv m 



2 + BC 2 + C 3 

S " v 

1PN 2PN 3PN 



(1 + B + C + BC + C + BC + C +0(4) • • •) , (2.9) 
dr r z ' ' ' • 



where the post-Newtonian order of the relevant terms is indicated by the PN label under the braces. These equations 
describe the truncated post- Newtonian model Q. Before proceeding ahead, we must verify that the above truncation 
is consistent and meaningful in the pres ent pr oble m of neutron stars. For a range of models, we evaluate the successive 
combinations that app ear in Eqs. ( pT^ ) and (2.9), using the exact TOV equations. Figure [l] shows the combinations 



appearing in Eq. ( |2.8| ) for two such models with an equation of state of intermediate stiffness, one which is very close 
to the maximum mass limit and another, a little further down the stable branch. From the Figure it is clear that 
numerically the successive combinations are perturbatively smaller to be consistently denoted as 1PN, 2PN and 3PN 
contributions. 

III. PADE APPROXIMANTS TO THE TRUNCATED POST-NEWTONIAN (TAYLOR) MODELS 

Given a Taylor series of order n in one expansion parameter x 
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S n (x) = l + a t x-\ ha„/, (3.1) 

its Pade approximant is the ratio of rational functions, 

= (3-2) 

satisfying the condition 

T n [P fe m (a;)] = S n (x); k + m = n. (3.3) 

In the above, N m and Dfc are polynomials in x of order m and A: respectively and T n denotes the operation of 
expanding a function in Taylor series upto order n. 

A convenient choice of Pade approximants is its diagonal (or nearly diagonal) continued fraction form |s| i.e., P™ 
when n = 2m (even) and P™+i when n — 2m + 1 (odd). For instance, 

Pl = > (3-4) 

1 H 



1 + c 3 x 

where the Cj 's are determined in terms of the Taylor coefficients by 

-, a\-a 2 aia 3 - a\ 

c = 1, ci = -ai, c 2 = , c 3 = — —J . (3.5) 

a\ a\[a\ - a 2 ) 

To construct the Pade approximant of the truncated post-Newtonian models discussed in the previous section, we 
must proceed carefully. In the problem under discussion, in general, there are three independent small parameters 
A, B, C, of the same order of smallness. Given this system of approximate differential equations containing three 
independent small parameters, can we construct an associated Pade like approximation with faster and improved 
convergence? Assume for a moment, that since A, B, C are of the same order of smallness, we decide to use as the 
variable of expansion the parameter B that takes the maximum value in the entire range of integration. The first 
order Taylor term A + B + C would be rewritten as B(l + A/ B + C / B). If we treat only B as the independent variable 
then its associated coefficient would be 1 + A/B + C/B. If A/B and C/B have only a weak dependence on r, we can 
indeed use a simpler Pade form with only one variable and this should suffice. We call this the one parameter Pade 
form for the TOV equations. In this case, the relevant equations at 2PN are the following:^] 



dm 
dr 

dp mpt 1 



47rr 2 pt , (3.6a) 

(3.6b) 

(3.6c) 



dr r 2 Ti 
dv ml 
dr r' 2 Q 2 ' 



A + B + C- AB-BC-CA-C 2 

Wher6 ' ^ 2 = A + B + C + A 2 +AB + CA + B 2 + BC ' (3 ' 6d) 
B + C-BC-C 2 

and < 52 = B + C + B 2 + BC ■ (3 - 6e) 

(The more involved equations^] at 3PN are listed in the Appendix). To examine more quantitatively, the validity of 
the above treatment, we consider, as before, a few models and compute the values of A, B and C and the associated 
ratios of A/B and C/B. The results are displayed in Figure From the figure it is clear, that as required, A, B, C 



1 Identical equations obtain if one uses any of A, B or C as the expansion parameter. It is also equivalent to introducing 
by hand a small parameter say e in terms of which we (Taylor or Pade) expand, treating as coefficients the associated A, B, C 
dependence. 

2 The explicit forms of -7-2,02 and the corresponding 3PN forms in the Appendix are exibited for a nalyt ical completeness, 
our numerical computation however, we directly substitute Eqs. (3.4) and (3.5) in Eqs. (3.6b) and (p.6q). 
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are approximately of the same order of smallness. However, if one looks at the ratio of C/B, one finds that there is a 
regime where the r-dependence is not very weak. Consequently, we are wary of using only the usual straightforward 
one parameter Pade form discussed above and proceed as follows. We generalise the usual construction for the one 
parameter Pade to the situation where there exist, not one, but three small parameters. The advantage of using a 
three parameter form is that its 'coefficients' are pure numbers with no explicit 'r-dependence'. By treating A, B, C 
as independent variables, we avoid the explicit issue of the 'weak r-dependence' of the associated coefficients^] 

To this end, we start with the most general second order post-Newtonian accurate polynomial in three small 
parameters A, B and C as 

T 2 = 1 + t A A + t B B + t c C + t AA A 2 + t BB B 2 + t cc C 2 

+t AB AB + t BC BC + t CA CA. (3.7) 

The associated Pade approximant in continued fraction form may be chosen as 

Pi=i, (3-8) 



where 



p 2 = 1 + ElA + 

1 + pnA + p 12 B + p 13 C 



P2B P3 ° (3.9) 



1 + p 2 iA + P22B + P23C 1 + p 3 i A + p 32 B + p 33 C 



By matching coefficients of the Taylor expansion of P/, Eq. (3^), and T2, Eq. (3/7), the three coefficients p\, P2 and 
p 3 are uniquely determined as, 

Pi = -t A ; P2 = -t B ; p 3 = -tc- (3.10) 

The remaining nine coefficients Pij , where i, j = 1,2, 3, are not uniquely determined; they are solutions to the following 
(under-determined) system of six equations 

(3.11a) 
(3.11b) 
(3.11c) 
(3. lid) 
(3. lie) 
(3.11f) 



t A (t A 


-Pll) 


— t AA , 




t B (t B 


-P22) 


= t BB , 




tc{tc 


-P33) 


= tec , 




t A P\2 + 


t B P2\ 


= 2t A t B 


— t AB , 


t B P2Z 4 


tcP32 


= 2t B tc 


~ tsc 1 


tcVzi 4 


t A Pl 3 


= 2t c t A 


- tcA ■ 



If t A ^ 0, t B ^ 0, tc ^ 0, pn, P22 and p 33 , can also be uniquely determined using Eqs. ( |3.11a| ), (|3.11h| ) and (|3.11c| ) 
respectively. We have 

Pn = t l^, (3.12a) 
t A 

t B ~ t BB 

P22 = , (3.12b) 

t B 

t 2 c - t C c /„ ir> , 

P33 = -^-7 ■ (3.12c) 

tc 



T he rem aining six non-diagonal terms cannot be uniquely determined from the remaining three equations Eq. (3. lid) 



(3.11i) without further input. The natural requirement that all Pade coefficients Pij contributing to a particular 



Taylor term, contribute equally, leads one to the following symmetry choice 



3 Indeed, one cannot relax this requirement in our three parameter construction either; however, this is only implicit in the 
fact that A, B, C are of the same order of smallness. 
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Pl2=P21, P23 = P32, Pl3=P31- 

One can then uniquely determine all the required coefficients and finally obtain 



(3.13) 



^A^B — ^AB 

P\2=P2\ = — : — — , if t A + t B ^0, 



P23 = P32 
P31 = Pl3 



t A ^ 


-tB 


2t B tc 


— tBC 


t B H 


-t c 


2t C t A 


- tcA 







, if t B +t c ^ 0, 
, if t c + t A ^0. 



(3.14a) 
(3.14b) 
(3.14c) 



Since the form of Eq. (2.8) is equivalent to Eq. (3.7) and Eq. (2.7) has a less general form containing only 2-parameters 
B and C, the second order Pade approximant to the second order truncated TOV equations may be written as 



dm o 

— — = 471T p t , 
dr 

dp mp t 1 
dr r 2 F2 

dv ml 



(3.15) 
(3.16) 
(3.17) 



where F2 is given by Eq. ( |3.9| ) and G2 is obtained by substituting A — in F2 . A comparison of Eqs. (2.8) and (3.7) 
yields 



t A — t B — tc — t AB — t B c — tc A — tec — 1. 

t BB = t AA = 0, 



so that the pij (Eqs. ( 3.12a )-( 3.12c ) and ( 3.14a )-( 3.14c )) in this case become 



Pll =P22 = 1, P33 = 0. 

1 

P12 = P23 = P13 = „ 



The associated forms of F2 and G2 are finally given by 

.4 



Fo = 1 — 



D 



C 



1 + A +M±c 1 + b+ a±c 



A+B 



G 2 



D 



C 



1 + B + 



(3.18a) 
(3.18b) 



(3.19a) 
(3.19b) 

(3.20) 
(3.21) 



Equations ( 3.20| ) and ( 3.21 ) define the three parameter 2PN Pade model associated with the 2PN truncated model 
discussed in S. The 3PN three parameter Pade model may be similarly obtained but since the expressions are 
lengthier, we list them in the Appendix together with the one parameter 3PN Pade model. 



IV. RESULTS AND DISCUSSION 

We consider the polytropic equation of state p — Kp T choosing the same models as studied by Shinkai i.e., 
r = 5/3, 2 and 3 with K (the polytropic constant) = 4.35, 10 2 and 10 5 , respectively. Before presenting the final results, 
we have compared the performance of the one parameter Pade models with three parameter Pade models and display 
in Figure [?] a typical comparison. From the Figure it is clear that both models yield similar numerical performance 
and the difference, if any, appears in the region where anyway these models are not recommended. Consequently, all 
our Figures refer to the one parameter Pade models though we have verified that the three parameter Pade models 
also yield similar results. 

In Fig. ^, we plot mass (in solar mass units) as a function of central density. The different curves represent the 
exact TOV, the 2PN and 3PN truncated (Taylor) models (henceforth 2-Taylor, 3-Taylor etc; for more details, see 
H) and our new Pade approximated models. We do not plot curves corresponding to 1PN truncated Taylor and the 
associated first order Pade model as they are very far away from the exact GR curve and thus evidently inadequate. 
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It is very clear from the figure that for all models to the left of the (first) maximum, which represents the stable 
branch of the mass-density curve, the Pade models do extremely well as compared to Taylor models. In this regime 
the third order Pade (now onwards 3-Pade) curve (dashed line) is very close to the exact GR curve (solid line). The 
2-Pade (long-dashed line) curve is not only better compared to the 2-Taylor one (dot-dashed) but even far better than 
the 3-Taylor curve (dotted line). In the unstable branch i.e., to the right to the maximum, the Pade approximation 
starts becoming bad. This feature is more severe when the stiffness of the considered equation of state and the order 
of the Pade approximants are higher^ (Fig. |^b and c). However, the stable branch is approximated well by Pade 
models even for a very stiff equation of state (r = 3, Fig. ^c), where the approximation breaks down near (2.75 x 
10 15 g/cm 3 ) the maximum mass limit (1.85 x 10 15 g/cm 3 ). Studies of mass-radius curves and v c (the exponent of 
the metric component gtt at the center of the star) vs central density curves also show a similar dependence on the 
stiffness of the chosen equation of state and on the central density. In Fig. |B| and Fig. we plot mass vs radius and 
v c vs central density curves respectively, for an equation of state with intermediate stiffness r = 2. The deviation of 
the 3-Pade solutions from the exact GR ones correspond to the configurations lying on the unstable branch. 

The studies of these equilibrium configurations imply clearly that the stable general relativistic TOV configurations 
are approximated very well by Pade models as compared to Taylor truncated models, even if they do not perform so 
well in the unstable branch. To study this further, we next numerically evolve these five initial configurations - Exact 
TOV, 2-Pade, 3-Pade, 2-Taylor and 3-Taylor truncated models - in time and qualitatively compare their evolution. 
We use a spherically symmetric general relativistic hydrodynamical code for this purpose jjj). The code uses polar 
slicing and radial gauge. The spacetime is evolved using the ADM formalism and the hydrodynamic evolution is based 
on a high resolution shock capturing |l(]]- [jll| scheme. The grid boundary is fixed at about four times the radius of 
the TOV model. We choose a model with central density 1.28 x 1CV 3 (geometrised units), r = 2, k = 10 2 and then 
evolve it upto 10 ms (though the code runs for a much longer time, 10 ms evolution is sufficient for our analysis). In 
Fig. ^ the evolution of central density using various initial data is displayed. We find that the 3-Pade curve follows 
the TOV evolution curve very closely. Even the 2-Pade evolution is much better than the 3-Taylor and the 2-Taylor 
evolutions. The oscillations in the central density are due to the numerical truncation errors which introduce non-zero 
radial velocity. These truncation errors act as the small perturbations on a stable spherically symmetric configuration 
and give rise to the radial pulsation modes of the system. A Fourier transform of the density or radial velocity 
evolution can be used to extract these pulsation modes Jl2]| and we hope to return to this in a subsequent work. In 
Fig. H we plot the metric components gtt and g rr over the grid for t — (initial time) and t = 10ms (final time). To 
compare the global performance of the initial data used (the exact TOV and the various truncated models), the norms 
of the Hamiltonian constraints (maximum value at a given time step) and l\ (average of the absolute value) are 
evaluated and shown in Fig. ||. These figures, once again, confirm the superior performance of Pade approximants over 
Taylor approximants of the same order. Though we do not display them, we get similar results for the other stable 
configurations both for relatively softer and stiffer equations of state. A final comment on the models in the unstable 
branch of the mass-density curve: the truncation errors are enough to trigger a collapse and make the system unstable 
on a time scale so short that a comparison of various approximations is not possible. A more detailed analysis would 
be needed in this regard. 

To conclude: Detailed studies of equilibrium configurations of single neutron stars and their evolution indicate, that 
in the stable branch, the second order Pade model converges to the exact general relativistic model even better than 
the straightforward third order truncated (Taylor) PN model. Both the simpler one parameter Pade form and a more 
involved three parameter Pade form exhibit similar improvement over the Taylor models. The Pade models are thus 
quite robust and controlled and perform better than the simpler Taylor truncated models. It is better to use initial 
data obtained from a Pade approximant to the Taylor model than initial data from a straightforward post-Newtonian 
truncated model of the same order. This feature should be generic and extend to binary neutron stars and black holes 
[especially since a useful simplification in a two-body problem is via a reduction to an equivalent one-body problem] 
and prove useful in numerical studies of such systems in the future. 



4 While computing the PN coefficients in Eq. (2.8) for the models which fall in the unstable branch, we notice that the values 
for the 1PN terms become even greater than one, indicating that neither the Taylor truncation nor the Pade expansion is 
reliable in this region. For the stable branch models these coefficients are always lesser than one as shown in Fig. [l] 
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APPENDIX A: 



The 3PN truncated model in general has the form 

T 3 



l + t A A + t B B + t c C + 
t AA A 2 + t BB B 2 + t cc C 2 + t AB AB + t BC BC + t CA CA + 
t AAA A 3 + t BBB B 3 + tcccC 3 + t ABC ABC + 

t AAB A 2 B + t BBC B 2 C + t CCA C 2 A + t ABB AB 2 + t BCC BC 2 + t AAC A 2 C . 



The associated 3PN Pade may be written as 



1 



where F3 = 1 + 



PiA 



I _|_ EllA _|_ P12-B _|_ P13C 
S11 S12 5i3 



+ 



P2B 



P21A _j_ P22-B _l P23C 



1 + 



S21 



S22 



+ 



S23 



1 + 



P3iA P32B j_ p 33 C ' 



S31 



S32 



+ 



(Al) 



(A2a) 



(A2b) 



•S33 



where Sij = 1 + pijiA + "pijiB + p^C ', i,j = 1 ... 3. The analysis discussed in Sec. Ill at 2PN may be repeated 
at 3PN after making the following natural symmetric choice implied by the requirement that all Pade coefficients 
contributing to a particular Taylor term contribute equally: 

(A3a) 
(A3b) 
(A3c) 

i = P312 = P321 , (A3d) 

(A3e) 
(A3f) 
(A3g) 

The remaining 10 independent 3PN Pade coefficients, uniquely determined by 1PN, 2PN and 3PN Taylor coefficients, 
are given by 



P112 


= P121 


= P2\\ , 


P122 


= P22\ 


— P2\2 , 


f>113 


= P131 


= P3U , 


P123 


= P132 


= P23\ = 


P313 


= P331 


= P133 , 


P322 


= P232 


= P223 , 


P332 


= P233 


= P323 • 



Pill = 



P222 = 



P333 = 



P112 = 



-t AA + t AAA t A 



t A (t A 2 


— tAA) 


—t BB 2 + 


t BBB t B 


t B (t B 2 


— t BB ) 


-tec 2 + tccctc 


tc (tc 2 


-tec) 



(A4a) 
(A4b) 
(A4c) 



-t A 4 t B + 2 t B 2 t A + (-t B 3 + (-2 t AB + 2 t AA ) t B + t AAB ) t A 
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+ (-2 tAA^AB + 2 tAAB^B + 2 t B 2 tAB) ^A - 2 tAA^B^ + t AAB^B^ 
—tAB{2tAA+tAB)t B | 

(t B + iyl) 2 (2 t A tB - tAA - tAB + tA 2 ) j , 
{-t B 2 - 2 t BB ) t/i 3 + (2 t B tAB + 2 t B 3 + tABB) tA 2 

(-is 4 + (2 tes - 2 tA B ) t B 2 + 2 tA BB t B - tA B it as + 2 t BB )) tA 

2 tA B t BB t B + tABBtB 2 | 

(ifi + ^a) 2 (tB 2 + 2 tAt-B - tAB - tBB) j , 

-tct4 4 + 2 t^c 2 + (-tc 3 + (2 t/tA - 2 t^c) ic + i/tAc) tA 2 
(-2 tyi^ylC + 2 t c 2 tAc + 2 t/tAc^c) *a — 2 tAAtc 3 + ^yiAC^C 2 
(*yic + 2t A A) tc 



{tA + tc) 2 (tA 2 + 2 t c t A - t A A - tAc) 

(-2 t B 2 t c + (-2 t c 2 - 2 t B c) t B -2 t BC tc) t A 3 
+ (-2 t c t B 3 + (2 t BC + 2t AC + 12 t c 2 ) t B 2 

+ (-2 t C 3 + (2 t A B -UbC + 2 t A c) t C + tABc) t B 
+ (2 tAB + 2 tsc) tc 2 + tABctc tA 2 



+ ((-2 tc 2 - 2 i^c) is 3 + (-2 t c 3 + (2 t^s + 2 t BC - 4 t AC ) t c + *abc) is 2 

+ ((2 tAC + 2 t BC - 4 t,4 B ) tc 2 + 2 tABctc - 2 tAB (t B C + *Ac)) *B 
-2 t C 3 tAB + tABC'tc 2 - 2 i/tC {tAB + tBc) tc^j tA 

-2 t c tB 3 tAc + ((2 i/tc + 2 i/ts) *c 2 + t ABC tc) t B 2 

+ (-2 t C 3 tAB + tABCtc 2 - 2 t BC {tAB + t A c) t C ) t B 



{tA + tc) {t B + t C ) {tB + tA) (2 tAtB + 2 t B t C + 2 t C tA - tAC - tAB - t B c) 

(-2 t C c - tc 2 ) tA 3 + (2 t c 3 + 2 t C t A c + tAcc) tA 2 
+ (^tc 4 + (2 t C c - 2 tAc) tc 2 + 2 t AC ctc - tAc {2tcc + t A c)) tA 

+tACctc 2 - 2tActcctc I 

^{tA+tcf (2t C t A -tcC-tAC+tc 2 )^ , 

I —tctB 4 + 2 t B 3 tc 2 + (—tc 3 + (2 tss — 2 isc) tc + t BB c) t B 2 
+ (—2 t BB t B c + 2 tc 2 t B c + 2 t BB ctc) t B 



P233 



— 2 tBBtc 3 + tBBC'tc 2 — tBC (2 t BB + t B c) tc | 

^{tB+tc) 2 (t B 2 + 2t B t c -t BC -t BB )^ , 

| ( _ 2 tec - tc 2 ) t B 3 + (t B cc + 2 t c 3 + 2 t BC tc) t B 2 

+ (— tc 4 + (2 tec — 2 t B c) tc 2 + 2 t B cctc — t B c (2 1 cc + isc)) is 



(A4i) 



+tscc^c —2t B ctcctc 



} 



(t B +t c ) 2 (2tijtc-icc+ic 2 -ti3c)| • (A4j) 

[The above solution obtains if none of the factors in the denominators of the above expressions are vanishing.] 
For the F 3 associated with the TOV equation at 3PN order, we have 

tcCC = tABC = tAcc = tscc = 1 j (A5a) 

tAAA = tBBB — tAAB = t-ABB = t B BC = t-AAC = . (A5b) 
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The solution for P333, Eq. (A4c), is not applicable, since tc — tec = 1- One finds that p 333 is indeterminate. We 
choose it to be zero. The function F 3 and G3 then reduce to 
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(A7a) 
(A7b) 

(A7c) 

(A7d) 



r -i B c 
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8 Z| > 1 8 4/ 



A> = 1 + 
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2(l-f "€) 



(A8a) 
(A8b) 
(A8c) 



9 



The corresponding one parameter Pade approximant at 3PN order on the other hand is given by 
F 3 = {A 3 + (1 + B + C)A 2 + [B 2 + (1 + C)B + C]A + B 3 + (1 + C)B 2 + BC}' 1 

[B 2 + (C-1)B-C+ l]A 2 + [{C - l)B 2 + (-2C + C 2 + 1)5 - C 2 + C]A + 



(1- C)B l + (C-C J )Bj, (A9a) 
(l-C)B 2 + (C-C 2 )B 

53 ~ B3 + (1 + C)B2 + BC ■ (A9b) 
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FIG. 1. Post-Newtonian coefficients appearing in Eq. ( |2.8| ) as a function of the radial distance (in units of km). The models 
are for a polytropic equation of state with intermediate stiffness (F = 2) with the central density p c = 8 x 10 14 g/cm 3 (panel 
(a)) and 1.85 x 10 15 g/cm 3 (panel (b)). The latter is a model which is very close to the maximum mass limit. 
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FIG. 2. Parameters A, £?, C and the ratios of parameters A/B,C/B as a function of radial distance. From the panels (c) 
and (d), we see the weak dependence of the ratio C '/B on r (the radial distance), though A/B is almost constant all through, 
except near the surface. 
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FIG. 3. A comparision of the two different Pade aprroximants constructed here: the one parameter Pade in panels [d-f] and 
three parameter Pade in panels [a-c] . It is clear that in the stable region both models give very similar results. 
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Central density (g/cm ) 

FIG. 4. Total mass (in units of Mq) as a function of central density (p c ) (in units of g/cm 3 ) for three values of T. Displayed 
here and in all subsequent figures are the one parameter Pade models. The three parameter Pade models also give very similiar 
results. Pade approximants do extremely well as compared to Taylor truncated models upto maximum mass limit i.e., for all 
the stable TOV models (models to the left of the first extrema). For the behaviour in the region beyond this extremum, see 
the discussion in section 4. 
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r = 2 





Radius (km) 



FIG. 5. Total mass (M© units) as a function of radius (km) for the intermediate value of F discussed in Fig. |^ i.e., r = 2. 
The deviation of the 3-Pade curve with respect to the GR curve for smaller radii models correspond to models (with high 
values of the central density) that fall in the unstable branch of the mass-central density curves. Pade models should not be 
used beyond the maximum mass limit. 
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FIG. 6. The exponent v of metric component g tt at r = as a function of central density. 
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FIG. 7. The evolution of central density for exact TOV, Pade models and Taylor truncated TOV initial configurations with 
the same value of central density p c = 1.28 x 1CP 3 and r = 2; k = 10 2 . The numerical truncation error triggers the dynamical 
evolution and shows the pulsation which corresponds to the physical pulsation modes of a stable spherically symmetric model. 
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FIG. 8. gtt and g rr over the grid at time t = Oms ((a) and (b), respectively). The insets zoom into particular portions of 
the grid (near the center for g u and near the surface for g rr ) to highlight the Pade behaviour more clearly, (c) and (d) are the 
same as (a) and (b), respectively, but at t = 10 ms. The model is the same as in Fig. ^. 
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